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O-lGeneral properties of accretion onto isolated stellar mass black holes in the Galaxy are discussed. An analysis of plasma internal 
^ fnergy growth during the infall is performed. Adiabatic heating of collisionless accretion flow due to magnetic adiabatic invariant 
j_i conservation is 25% more efficient than in the standard non-magnetized gas case. It is shown that magnetic field line reconnections 
in discrete current sheets lead to significant nonthermal electron component formation, which leads to a formation of a hard (UV, 
d }f-ray, up to gamma) , highly variable spectral component in addition to the standard synchrotron optical component first derived 
1 'by Shvartsman generated by thermal electrons in the magnetic field of the accretion flow. Properties of accretion flow emission 
j variability are discussed. Observation results of two single black hole candidates - gravitational lens MACHO- 1999-BLG-22 and 
radio-loud x-ray source with featureless optical spectrum J1942+10 - in optical band with high temporal resolution are presented 
and interpreted in the framework of the proposed model. 

m : 
in ■ 

Key words: Black holes, Infall and accretion, Photometric, polarimetric, and spectroscopic instrumentation 
. PACS: 97.60.Lf, 98.35.Mp, 95.55. Qf 
OS ' 



O 

o 



1. Accretion onto isolated stellar- mass black holes 



K^J Even though more than 60 years have passed since the 
theor etical prediction of black hole s as an astrophysical ob- 
d jects ( OppenheimerfcSnvder , 1939h in some sense they have 
not been discovered yet. 

The features of the black hole - the compactness (the size 
for the mass M is close to the Schwarzschild one r g = 2G c ^ 1 
and the mass larger than 3Mq - are the necessary, but not 
the sufficient ones. The deciding property of the black hole 
is the presence of the event horizon, instead of with the 
usual surface. To decide whether it is present in a given 
compact and massive object it is necessary to detect and 
study the emission generated very near the event horizon. 

It is a very complicated task that cannot be easily per- 
formed in x-ray binaries and AGNs due to high accretion 
rate and so the high optical depth of the accreting gas. At 
the same time, single stellar-mass black holes, which ac- 
crete interstellar medium of low density (10 -2 — lcm~ 3 ), 
are t he ideal case for de tection and study of the event hori- 
zon ( Shvartsmanl . 1971 ). 

The analysis of existing data on possible black hole 
masses and velocities in comparison with the interstellar 



medium structure shows that in the majority of cases in 
the Galaxy (> 90%) the accretion rate rh = M (?/L e _dd 
can't exceed 10" 6 - 10" 7 (|Beskin fc KarpovLl2005t ). 

Black holes usually move supersonically (at Mach num- 
bers 2-3). And only in cold clouds of the interstellar hydro- 
gen (n ~ 10 2 — 10 5 , T ~ 10 2 ) and at low velocities of motion 
(< 10 km/s) the accretion rates may be high and black hole 
luminosities may reach 10 38 — 10 40 erg/c. For typical inter- 
stellar medium inhomogeneity the captured specific angu- 
lar momentum is smaller than that on the BH last stable 
orbit, and the accretion is always (with the exception of the 
case of a bl ack hole in a dense molecular cloud) spherica lly- 
symmetric (jMiifcTotanil . 120051: iBeskin fc Karpovl l2006t) . 

The accreting plasma is initially collisionless, and it re- 
mains so until the event horizon. The electron-electron and 
electron-ion free path A ~ 2.4 • 10 3 T 2 n~ 1 even at the cap- 
ture radius is as high as ~ 10 12 cm. Only the magnetic 
fields trapped in plasma (the proton Larmor radius at r g 
is 10 cm) make it possible to consider the problem as a 
quasi-hydrodynamical one; it is only due to the magnetic 
field that the particle's momentum is not conserved, allow- 
ing particles to fall towards the black hole. In addition, the 
magnetic field effectively "traps" particles in a "box" of 
variable size, which allows us to consider its adiabatic heat- 



ing during the fall; a correct treatment of such a process 
shows that for magnetized plasma such heating is 25% more 
effective than for ideal gas due to the conservation of mag- 
netic adiabatic invariant I = , whe re p t is the tangen- 



tial c omponent of electron momentum dBeskin fc Karpov , 
2005h . Therefore, the plasma temperature in the accre- 
tion flow grows much faster and electrons become relativis- 
tic earlier - R re i sa 6000 in contrast to R re i ~ 1300 in 
Bisnova tvi-Kogan fc Ruzmaikinl (|1974j ) and R re i w 200 in 



Ipser fc Pried ( 1982t ). The accretion flow is much hotter, 
and our estimation of " thermal" luminosity 

L = r]Mc 2 = 9.6 ■ 10 33 M^ n?(y 2 + c 2 )^ 3 erg/s (1) 



is significantly higher than those of Ipser fc Pried ( 19821 ) 



Lip = 1.6 • 10 32 M 3 n 2 (F 2 + c 2 s )^ erg/s (2) 



and lBisnovatvi-Kogan fe Ruzmaikinl (|1974l ) 

L BK R - 2 ■ 10 33 M 3 n 2 (y 2 + c 2 )^ 3 erg/s, 



(3) 



while the optical spectral shape is nearly the same. The ef- 
ficiency of accretion as a function of accretion rate is shown 
in Fig. [5] 

The basis of our analysis is the assumption of the energy 
equi partition in the ac cretion flow (Shvartsman' theorem, 



sec; 



Shvartsman (Il97lh . The straight consequence of this 



assumption is the necessity of exceeding magnetic energy 
dissipation with the rate 



dE 



_ 4 »B 2 hvB 2 _3vB 2 
dVdt~ r 8?r ~~ 2 r 8tt ~ 2 r 8tt ' 



(4) 



In the pre vious accretion models the d i ssipati on goes 
continuously) Bisnovatvi-Kogan fc Lovelacd . 1997? ) in the 
turbulent flow and its mechanism is not examined in de- 
tails. We considered conversion of the magnetic energy 
in the discrete turbulent current sheets (|Pustilnikl . fl997f l 
as a mechanism providing such dissipation. At the same 
time, various modes of plasma oscillations are generated 
(ion-acoustic and Lengmur plasmons mostly), magnetic 
field lines reconnect and are ejected with plasma from the 
current sheet and electrons are accelerated. The latter ef- 
fect is very important for the observational appearances 
of the whole accretion flow. The beams of the accelerated 
electrons emit its energy due to motion in the magnetic 
field and generate an additional (in respect to synchrotron 
emission of thermal particles) nonthermal component. 

Current sheet spatial scales are usually much smaller 
than the whole accretion flow one and, therefore, the frac- 
tion of the accelerated particles is small, and so the to- 
tal (significantly nonthermal) electron distribution f(R, 7) 
may be considered as a superposition of the purely thermal 
one ft(R, 7) for the background flow particles and purely 
nonthermal f n t{R,l) for the ones accelerated in the cur- 
ren t sheets (this i s the approach known as " hybrid plasma" 



see 



Coppil (|l999f) 1. Also we assume that for low accretion 



rates the non-adiabatic heating and radiative energy losses 



do not change the shape of the thermal component distri- 
bution (while changing its mean energy), then 



f(R,j) = ft(R,~/)+(fnt(R,l). 



(5) 



Note that this distribution is not normalized, and only its 
shape has physical meaning. So, for example, the ratio of 
nonthermal to thermal electron densities at some radius R 
may be expressed as 



n nt {R) = CfntjR) 
n t (R) ft(R) 



(6) 



where f n t (R) and ft (R) are integrals of the corresponding 
distribution functions over the range of 7. 

Thermal particle distribution function may be written as 



Vr a 



(-?)• 



(7) 



where the usual dimensionless expression for temperature 
r = kT/m e c 2 is used. This gives a Maxwellian local energy 
distribution and radial density slope p ex R~ 3 ^ 2 . 

To get the temperature profile of the accretion flow let's 
take into account that the fraction £ ~ 0.1 of the magnetic 
field energy during its dissipation is converted to the energy 
of the accelerated particles. Then, according to (U]), the 
non-adiabatic heating rate may be expressed as 



* = (i-0— = (i-0~- 

V ^'dVdt V W 2r87r 



(8) 



The main mechanism of radiative losses at low accretion 
rates is synchrotron radiat ion. Its rate may be written as 
(|LightmanfcRvbickil . fl979h 



A sync = gCr T c7"— n. 



8tt 



For a Maxwellian distribution 
12 [ _^_V = 12r 2 . 



(9) 



(10) 



For relatively large accretion rates (m > 10 ~ 6 ) it is nec- 
essary to take into account the cooling due to inverse Comp- 
ton effect, which converts the part of hot electrons energy 
to the hard emission. 

The rate of inverse C ompton energy losses is re lated to 
the synchrotron one as (|LightmanfcRvbickil [l979h 



A-compton — A-sync 



(ii) 



B 2 /8tt 

Taking into account non-adiabatic terms, the energy 
balance equation per particle may be written for a non- 
relativistic region of the flow (R > R re i ) as 

d e _ 5 e 1 A - $ 
dr n Anr 



v n 



(12) 



For a non-relativistic electron gas e = |nfcT, so 
dr_ = _5t 

dR 4R ?j 2 m e 



(13) 



Fig. 1. Decomposition of a single black hole (with the mass IOMq) 
emission spectrum into thermal and nonthermal parts. The accretion 
rate is 1.4 ■ 10 10 g/s, which corresponds to fa = 10~ s . 



while for the relativistic one (e — 3p — 3nkT) 



— = -t- - (i - ^ — IHpr-Z 
dR 8R [ 4 m e 



4 uip rhr 2 
3 m, R 2 



(14) 



where a 2 — -| (the equipartition condition). 

The boundary condition is r(R re i) = 1. 

An analytical solution of this equation in general is very 
difficult, but for low accretion rates we may neglect the 
influence of radiative losses and get a solution in the form 



t(R) = (1-0 2 4^K- 1 
3 m e 



2a 2 m„ 

— - 

SKrel m e 



This expression may be substituted into (|7|) to get the final 
expression for the thermal electron distribution. 

The energy evolution of single nonthermal electron is 
described by 



eh, 
dR 



1 tn-p . 7 2 5 7 

a^i 771 ^ 2 " ~ &r' 



(16) 



where the first term corresponds to synchrotron losses and 
the second one to adiabatic heating. 

Analysis of the e volution of electron beam s ejected from 
the current sheets (jBeskin fc Karpovl 120051 ) gives the pos- 
sibility to compute the £ and f n t{R, 7), i.e. the fraction of 
thermal and nonthermal component contributions to the 
total one. 

The spectra of various components of the accretion flow 
emission have been computed taking into account rela- 
tivistic effects of emission propagation ne ar the black hole 
(jLightmanfcRvbickil . Il979llshapird . 1 19731) . They are shown 
in Fig. [5] for various accretion rates. In Fig. [3] the integral 
dependence of the thermal emission of accretion flow on the 
distance to the black hole is presented. The large fraction 
of it is generated inside the 2r g sphere, and so carries the 
information very near the event horizon. 

An important property of the nonthermal emission is its 
flaring nature - the electron ejection process is discrete, and 
typical light curve of single beam is shown in Fig. 21 The 
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Fig. 2. Spectra of the accretion flow onto the 10 Mq black hole for 
the various accretion rates. 




Fig. 3. Fraction of thermal synchrotron emission that comes from 
inside a given radius R 



light curve of each such flare has a stage of fast intensity 
increase and sharp cut-off, its shape reflects the properties 
of the magnetic field and space-time near the horizon. 

The black hole at a 100 pc distance (a sphere with 
this radius must contain seve ral tens of such objects, see 
( Agol fc Kamionkowskll2002l )) looks like a 15-25 m optical 
object (due to the "thermal" spectral component) with a 
strongly variable companion in high-energy spectral bands 
("nonthermal" component). The hard emission consists 
of flares, the majority of which are generated inside a 5r g 
distance from the BH. These events have the durations 



~ r g /c (~ 10~ 4 s), a rate of 10 3 -10 4 flares per second, and 
an amplitude of 2%-6%. The BH variable X-ray emission 
may be detected by modern space-borne telescopes. 

Optical emission consists of both a quasistationary 
"thermal" part and a low- frequency tail of nonthermal 
flaring emission. The rate and duration of optical flares 
are the same as X-ray ones, while their amplitudes are 
significantly smaller. Indeed, the contribution of nonther- 
mal component to the optical emission is approximately 
2 • 10~ 2 for rh = 10~ 8 — 10 -6 , so the mean amplitudes of 
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Fig. 4. Internal structure of a flare as a reflection of the electron cloud 
evolution. The prevailing physical mechanisms defining the observed 
emission are denoted and typical durations of the stages are shown. 
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Fig. 5. Efficiencies of the synchrotron emission of thermal and 
non-thermal electron components of the accretion flow. 

optical flares are 0.04%-0.12%, while the peak ones may be 
1.5-2 times higher and reach 0.2%. Certainly, it is nearly 
impossible to detect such single flares, but their collective 
power reaches 18-24 m and thus may be detected in obser- 
vations with high time resolution (< 10~ 4 s) by the large 
optical telescopes (jBeskin fc Karp"ov , 2005f ). 

Of course, the variability of the BH emission is related 
not only to the electron acceleration processes described 
here. Additional variability may be result of the plasma 
oscillations of different kinds, or other types of instabilities. 
The time scale of such variability may be from r g /c till r c /c, 

1. e. from microseconds till years. 

2. Search for isolated stellar mass black holes in 
optical range 

As it has been already noted, the main indirect ob- 
servational sign of the isolated stellar mass black hole is 
the absence of any spectral features in the emission of 
accreting pla s ma (e xcept for interstellar absorption ones) 
( Shvartsmanl Il97lf h Among such objects are the white 
dwarfs with featureless spectra (DC-dwarfs) and known 
localization (as they are usually selected as an objects 



with known proper motions) (iMcCookfeSionl . Il999h . On 
the other hand, in the last years by means of the cross- 
correlation of surveys of various wave bands (from radio to 
gamma) and the follow-up spectroscopic observations, the 
population of objects with featureless optical s pectra and 
without known localiz a tions, has been detected (jPustimikl . 



1977t iTsarevskv et all . 120051) . The latter include the ra- 



dio objects with continuous optical spectra (ROCOSes) 



Pustimikl.ll977h. Fo r them , as well as for some DC-dwarfs 



( Be skinfcMitronoval 1991), the upper limits on the fast 
optical variability as an observational manifestation of iso- 
lated stellar mass black holes ha ve been derived for the first 
time ( Shvartsman et al . 1989af bt). This idea became the 
basis of the observational programme for search for BH on 
Russian 6-m telescope - MANIA (Mu ltichannel Analysis 
of Nanosecond Int ensity Alterations) (jShvartsmanl . 11977 : 
iBeskinet aj|l997l) . 

Recently, some evidences appeared that the isolated stel- 
lar mass blac k holes may be among the unidentified gamma - 
ray sources (|Gehrels et a 1 120001 : La Palombara et al, 
The most promising cases are such where the indepen- 
dent estimation of the object mass is possible. One of such 
cases is the black hole in a binary system with a white 
dwarf. Such a binary may be detected by means of its pe- 
riodic brightness amplification on the t ens of seconds time 
scale due to gravitational microlensing (jBeskin fc Tuntsovl . 
l2002h . By using the technical equipment of SDSS (23 m limit 
in 6 square degree field) roughly 15 such objects may be 
detected in 5 years. It is clear that in such system it is easy 
to estimate the mass of the black hole, while there is no 
accretion from the white dwarf onto it, and the black hole 
interacts with the interstellar medium only. Another case 
when we may estimate the mass o f the black hole are the 
long MACHO microlensing events (jBennett et all . [2002) . 

Below we present the results of observations of two can- 
didates to the isolated stellar mass black holes. 



3. Observations of objects-candidates to isolated 
stellar-mass black holes at 6-m telescope 

Sample observations of the longes t micr olensing event 
MACHO 1999-BLG-22 (|Bennett et all l2002h . a most solid 
single stellar-mass black hole candidate, have been per- 
formed at the Special Astrophysical Observatory of 
RAS in 2 003-2006 in the f ramework of the MANIA 
experiment (jBeskin et al 1997T ). 

We used in observations the multichannel panoramic 
spectro-polarimeter (MPSP) based on posit ion-sensitive 
detector (PSD) with 1 11s time resolution (jPebur et all . 
l2003HPlokhotnichenko e"tall2003l) (see Fig.©. Such detec- 
tors use the set of multichannel plates (MCP) for electron 
multiplication, and multi-electrode collector to determine 
its position. PSD used in our observations has the following 
parameters: quantum efficiency of 10% in the 3700-7500 
A range, (S20 photocathode) , MCP stack gain of 10 6 , 
spatial resolution of 70 /im (0.21" for the 6-m telescope), 
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Fig. 6. Optical scheme of the multichannel panoramic spectro-po- 
larimeter (MPSP). 




Fig. 7. Multichannel panoramic spectro-polarimeter. 

700 ns time resolution, 7 • 10 4 pixels with 22 mm working 
diameter, and the 200-500 counts/s detector noise. The 
acquisition system used is the "Quantochron 4-480" spec- 
tal time-code convertor with 30 ns time resolution and 10 6 
counts/s maximal count rate. 

We combined our data with the archival observations 
of this region with XMM-Newton and Hubble ACS. We 
acquired the upper limits on the black hole x-ray emission 
on the level of 10~ 14 erg/s/cm 2 . 

The Hubble ACS optic al data allow to de-b lend the orig- 
inal source detected by (jBennett et a 1, 2002) on the level 
of / = 19. l m . We performed the /-band photometry of the 
three stars of the blend and placed an upper limits of / = 
19.5 m , / = 20.0 m and / = 21. 8 m on the /-band emission of 
the BH for the three possible cases. 

The parameters derived from the microlcnsing light curve 
don't allow to get the unique value for the distance (and 
so, the mass) to the black hole. Instead, it points to the 
two possible configurations of the microlensing system, one 
with the lens at 500 pc (with M — 136Af Q ) and the second 
with 2 kpc distance (M — 3.5Af ). We also studied the 



Fig. 8. MACHO-1999-BLG-22, observed with panoramic photometer 
with high temporal resolution of Russian 6-m telescope. Overplotted 
is the Hubble ACS image with the blend resolved. Green circle is 
the position of the lens error box. 
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Fig. 9. The prediction for the observational parameters of the MA- 
CHO-1999-BLG-22. The x-ray flux is in the 0.2-10 keV band. Over- 
plotted are the limits from the XMM-Newton and Hubble ACS 
archival observations. 

intermediate case (2 kpc distance, M = 27M©). 

The predicted observational parameters for these cases, 
computed according to our model, are shown in Fig.[9l with 
the observational limits overplotted. 

Our results rule out the case of the near and massive 
black hole. 

The analysis of data from 6-m telescope do not reveal any 
variable source. The upper limit on the variable emission 
component is placed on the B = 22 m level for the 10~ 6 - 
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Fig. 10. The optical Fourier power spectrum for the radio-loud x-ray 
source with featureless optical spectrum J1942+10. The upper limit 
for the variable emission component is 10% in 10 -5 -Is range. 

100 s time scale. 

The J1942+10 is the radio-loud x-ray object with fea- 
tureless optical spectrum with B ~ 18 m . We observed it 
for 40 minutes in photometric mode and placed the upper 
limits on the B-band variability on the 10% level over the 
10~ 5 -Is time scale (see Fig. ITO)) . 
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